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We discuss a many-electron Hamiltonian with Hubbard-like repulsive interaction and linear coupling 
to the phonon branches, having the Cu-O plane of the superconducting cuprates as a paradigm. A 
canonical transformation extracts an effective two-body problem from the many-body theory. As a 
prototype system we study the CUO4 cluster, which yields electronic pairing in the Hubbard model; 
moreover, a standard treatment of the Jahn- Teller effect predicts distortions that destroy electronic 
pairing. Remarkably, calculations that keep all the electronic spectrum into account show that 
vibrations are likely to be synergic with electronic pairing, if the coupling to half-breathing modes 
predominates, as experiments suggest. 
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The role of electron-phonon (EP) interactions in determin- 
ing the superconducting correlations in the Cu-O planes of 
cuprates is a very controversial issue. Possibly, the pair- 
ing mechanism has a predominantly electronic origin [1], but 
many high-Tc compounds exhibit a quite noticeable doping 
dependent isotope effect [2], suggesting that EP interactions 
could be important and should be included in the theory. 
In particular there is experimental evidence [3] that the half 
breathing Cu-O bond stretching mode at k — (n,0), (0, n) is 
significantly coupled with the doped holes in the supercon- 
ducting regime and its contribution may be relevant for the 
d x 2_ y 2 pairing [4] [5] [6]. 

The simplest way to include both strong electronic corre- 
lations and EP interactions is the Hubbard-Holstein model, 
where electrons are coupled to a local Einstein phonon. Much 
is known [7] [8] about the possibility of a superconducting 
phase in this model; further evidence for pairing can be ex- 
tracted by exact numerical diagonalization of small cluster 
Hamiltonians, calculating the effective pairing interaction de- 
fined as A(N + 2) = E(N + 2) + E(N) - 2E(N + 1), where 
N is the number of fermions in the system. A negative A is 
interpreted as an effective attraction and suggests the pres- 
ence af a bound pair in the ground state. Petrov and Egami 
[9] found A < in a doped 8-site Hubbard-Holstein ring and 
showed that the effect disappears once phonons are turned 
off. Conversely Mazumdar and coworkers [10] suggested that 
pairing in Hubbard clusters reported by some authors [If] 
was of doubtful physical interpretation due to the neglect of 
the lattice degrees of freedom and the Jahn- Teller (JT) effect. 
Consider even N, such that the N + 2-particle cluster hosts a 
singlet pair and the N + 1-particle ground state is degenerate. 
They argued that JT distorsion might cause a larger energy 
gain of the system with N + 1 particles, thus reversing the 
sign of A. This led to the conjecture that any A < due 
to an electronic mechanism is just a finite size effect, which 
vanishes for large systems like the JT effect does. 

Here we further investigate this issue by addressing the 
question if the recently proposed W — pairing avaliable 



in the Hubbard model [12] [13] [14] survives when the lat- 
tice degrees of freedom are switched on. In previous works it 
was shown that a class of Ca v -symmetric clusters exhibit the 
A < property, due to the binding of the so called W — 
pairs. This depends crucially on the symmetry and on the 
availability of degenerate states. Fermions in strongly cor- 
related systems can avoid the on-site repulsion and possibly 
form singlet bound states if there is enough symmetry; the 
fillings and the symmetry channels where the W = pairing 
can occur are determined in full generality by the W = the- 
orem [15]; these symmetries achieve the same result as high 
angular momentum and parallel spins in the Kohn-Luttinger 
[16] continuum approach. 

When lattice effects are introduced in this scenario, sev- 
eral questions arise. In the conventional mechanism, phonons 
overscreen the electron repulsion; what happens if electronic 
screening already leads to pairing? It is not obvious that 
the phonons will reinforce the attraction while preserving the 
symmetry. More generally, some vibrations could be pairing 
and others pair-breaking. To address these problems we use 
an extension of the Hubbard model in which bond stretch- 
ings dictate the couplings to the normal modes of the Ca v - 
symmetric configuration, generating a long-range (Frohlich) 
EP interaction. This is physically more detailed than the 
Hubbard-Holstein model, and does not restrict to on-site EP 
coupligs that would be impaired by a strong Hubbard re- 
pulsion. A further important drawback of the Holstein EP 
interaction in this context is an unphysically large polaron 
(and bipolaron) mass and low polaron mobility. This prob- 
lem is avoided using Frohlich-like phonons for modeling the 
Cu-O planes. It was pointed out recently [17] [18] [19] that 
even in the strong EP coupling regime the a long-range EP 
interaction removes the problem of polaron self-trapping. 

We start from the Hubbard model with on-site interac- 
tion U and expand the hopping integrals ti,j(ri, r 3 -) in powers 
of the displacements pi around a Ci v -symmetric equilibrium 
configuration 
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where a — x,y. Below, we write down the pf in terms of 
the normal modes q vv : pf = Y] S% v {i) q V v, where rj is the 
label of an irreducible representation (irrep) of the symmetry 
group of the undistorted system and v is a phonon branch. 

Thus, treating the Cu atoms as fixed, for simplicity, one 
can justify an electron-lattice Hamiltonian: 



H () + V tc 



(2) 



Here Ho = Hq + Hq is given by 

H = ^2 fi u V,v b v," b V,v + ^ t°j (l-j , Tj ) {c\ a Cj a + h.c), (3) 

where uj v ,v are normal mode frequencies. Moreover, let M 



denote the O mass, = \ v 



with A,, v numbers 



of order unity that modulate the EP coupling strength. Then, 
Vtot = V + W reads 
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the H Vt „ operators are given by 
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Consider the Hubbard model Hh = Hg + W involving the 
purely electronic terms of H e i-i att . We recall that W = 
pairs are defined as two-body singlet eigenstates of the Hh 
free of double occupancy; in the many-body problem they in- 
teract indirectly and can get bound [12] [20]. Here we wish 
to derive an effective interaction between the particles in the 
pair suitable for H e i-i a tt, by generalizing the canonical trans- 
formation approach of Ref. [20] . 

Denote the phonon vacuum by |0}} and the iV-particle non- 
interacting Fermi sphere by |$o(A r )); if we add a W = 
pair to ]$ (A i ')) ® |0)), the two extra particles, by definition, 
cannot interact directly (in first-order). Hence their effective 
interaction comes out from virtual electron-hole (e-h) excita- 
tion and/or phonon exchange and in principle can be attrac- 
tive. To expand the interacting (A r + 2)-fermions ground state 
1*0(^ + 2)), we build a complete set S of configurations in 
the subspace with vanishing z spin component, considering 
the vacuum state |$o(A0) <8> |0}} and the set of excitations 
over it. 

We start by creating W — pairs of fermions over 
|$o(A0)®|0)); we denote with |m)®|0)) these states. At weak 
coupling, we may truncate the Hilbert space to the simplest 
excitations, i.e., to states involving 1 e-h pair or 1 phonon cre- 
ated over the |m)<8»|0)) states. We define the |m)<g>|g}} states, 
obtained by creating a phonon denoted by q = (77, v) over the 
|m) ® |0}} states. Finally we introduce the \a) ® |0)) states, 



obtained from the \m) ® |0)) states by creating 1 electron-hole 
(e-h) pair. We now expand the interacting ground state in the 
truncated Hilbert space: 

|*o(JV + 2)) = J2a m \m) ® |0» + 

m 

+^a m , 9 |m) ® \q)) + ^a Q |a) ® |0)) (6) 

m,q a 

and set up the Schrodinger equation with energy eigenvalue E. 
A canonical transformation [21] decouples the higher ampli- 
tudes a m>9 and a a , and the a m become expansion coefficients 
over the W — pairs of the wave function of the dressed pair 
\ip). This obeys the Cooper-like equation H pa i T \(p) = E\tp) 
with an effective two-body Hamiltonian H v ^ r : 



H pair = {H + W + S[E]) 



(7) 



where S is the E— dependent effective scattering operator; the 
pairing problem must be solved self-consistently. It turns out 
that the effective scattering operator is 

or \ ^ W m , a W atm i \ - V m,m" V rL" ,m' 

b[h\ m , m , - E , _ E + 2^~e> — 1 W 

at a rr 

a m" ,q 



Here the primed quantites in the denominators are the 
eigenenergies of the \m) ® [0)), \m) ® \q)}, \a) ® |0)) states; 
they get renormalized in the decoupling procedure. The pair- 
ing criterion involves the properly renormalized Fermi energy 
£p (see ref. [13]); if the lowest energy eigenvalue E is such that 
E — 2ep + A with negative A, the dressed W = pair gets 
bound in the many-body interacting problem and the system 
undergoes a Cooper instability. 

As an illustrative application of the above pairing scheme, 
in this preliminary work we focus on CUO4, the smallest clus- 
ter yielding W — pairing in the Hubbard model. This re- 
quires 4 holes, (total number, not referred to half filling); such 
a doping is somewhat unrealistic, but larger C-iu-symmetric 
clusters and the full Cu02 plane also show W — pairing in 
the doping regime relevant for cuprates [12] [20]. Remarkably, 
in the pure Hubbard model, one can verify that A = A(4) at 
least at weak coupling [12], which demonstrates that A has 
the physical meaning of an effective interaction. CuCU rep- 
resents a good test of the interplay between electronic and 
phononic pairing mechanisms since we can compare exact di- 
agonalization results with the analytic approximations of the 
canonical transformation. CuCU allows only the coupling to 
phonons at the centre or at the edge of the Brillouin Zone; 
however, phonons near the edge are precisely those most in- 
volved [3] [4]. 
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FIG. 1. Pictorial representation of the ionic displacements in 
the eight normal modes of the Cu04 cluster, labelled according to 
the irreps of the Cav Group. 



Even in this small system the virtually exact diagonaliza- 
tions are already hard and the next C4„-symmetric exam- 
ple, the CU5O4 cluster [22], is much more demanding for the 
number of vibrations and the size of the electronic Hilbert 
space. The pure Hubbard CUO4 cluster with 0-0 hopping 
t x = yields [12] A(4) < 0, due to a couple of degen- 
erate W — bound pairs, in the Ai and B2 irreps of the 
Civ group; therefore in Eq.(8) we set the m — m' labels ac- 
cordingly. Starting with the C*4i>-symmetric arrangement, any 
displacement of the Oxygens in the plane can be analised in 
irreps, Ai, A2, Bi, B2, Ei, E2, see Fig.l. The EP coupling is 
expressed through two parameters dt and dt ox , defined as the 
derivatives of t and t ox with respect to the Cu-0 and 0-0 
bond lengths. The Cu04 cluster is also interesting as a test 
of the conjecture of Ref. [10]; we checked [21] that a standard 
JT calculation in which the degenerate ground state wave 
functions interact with the vibrations indeed predicts distor- 
tions that already at moderate EP coupling destroy W = 
pairing; here we wish to go beyond this approximation. 

According to Eq.(8), the phonon-mediated interaction, in 
the A\ channel is: 
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Thus, the two W = pairs behave differently: the B2 
binding energy is enhanced by phonons, while in the A\ sector 
the overall sign depends on A^ and lo v . It turns out that Ai 
and B2 modes are synergic to the W = pairing, while both 
longitudinal and transverse E modes are pair-breaking. The 
half-breathing modes that are deemed most important [4] [5] 
are Ai ± B\ combinations, but B\ is scarcely relevant. For 
the sake of argument, in the explicit calculations we took all 
the normal modes with the same energy erj = 
and A,, = 1. This sets the length scale of lattice effects £0 = 

^2M^ - 10_1 A where we used M = 2.7 x 10~ 26 Kg for 
Oxygens. 




FIG. 2. Analytical results of the canonical transformation: 
pair binding energy in the A\ and B2 sectors as a function of dt ox . 

Here we used = 1 for every mode, t = leV, t ox = 0, U = leV; 

0-1 

dtox is in units of eo/£o = leVx A , A is in eV. 



With this choice, pairing is enhanced in the Ai sector as 
well, albeit less than in B2, see Fig. 2, where the pairing en- 
ergy without phonons is A ~ — 20meV for both the W = 
pairs. In other terms, the vibrations split the degeneracy 
of the electronic ground state, effectively lowering the sym- 
metry like a nonvanishing t ox . Here again we stress that 
A = A(4) for both symmetries, generalizing the Hubbard 
model result. To go beyond the weak coupling case, we diag- 
onalized H e i-i at t numerically in the truncated Hilbert space 
with up to 5 modes at a time having vibrational quantum 
numbers n v < 3. In Figure 3a we included all vibrations 
except E\\ with increasing the EP coupling, the pair bind- 
ing energy steadily increases, in agreement with the trend in 
Fig.2; the unbinding mode E2 is dominated by the pairing 
ones. Note that the binding energy increases with dt ox as 
predicted by the weak coupling theory, but also with \dt\, 
suggesting that at strong coupling the Cu-O bond stretching 
grows in importance. 
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FIG. 3. A(4) as a function of dt ox for different values of dt. 

X v = 1 all the vibrations, except: A,, = for n = E\ (a); A,, = 

for r\ = E2 (b). n v =3 for each active mode; dt = —0.2 (diamonds); 

dt = —0.5 (triangles); dt = —1 (crosses). Here we used t = leV, 

0-1 

tox = 0, U = leV, dt ox and dt are in units of £o/£o = leVx A , 
A is in eV. 



In Figure 3b we included all vibrations except E^,. The 
E\ vibrations are strongly pair breaking, and A(4) is much 
smaller in magnitude than in Figure 3a; it becomes positive 
at intermediate coupling. It is likely that the couplings to the 
pair-breaking E modes are somewhat underestimated by the 
choice of A parameters in Fig. 3a and overestimated in Fig. 3b. 
Interestingly, however, at stronger coupling the pairing modes 
prevail again driving the system to pairing (bipolaron forma- 
tion). This resembles the trend found in Ref [9], although 
there a Hubbard-Holstein framework is used and since the 
system is one-dimensional the electronic pairing does not oc- 
cur. 

These results yield several indications. First, the catas- 
trophic distortions predicted by the JT approximation are not 
borne out by more complete approaches involving a broader 
spectrum of electronic states. This interesting result arises 
since the electronic ground state belongs to a narrow mul- 
tiplet; the standard treatment of the JT effect singles out 
the lowest level, which is justified only at weak EP coupling. 
Pairing prevails also at strong coupling in part of the param- 



eter space, in the symmetry channels where W = pairs 
occur. The correct trend is predicted by the canonical trans- 
formation approach, which also explains the pairing or pair- 
breaking character of the modes. In particular it is found that 
the half-breathing modes give a synergic contribution to the 
purely electronic pairing; since they are believed to be mainly 
involved in optimally doped cuprates, our findings suggest 
a joint mechanism, with the Hubbard model that captures a 
crucial part of the physics. This scenario was also drawn in the 
context of an extended t-J model, where the half-breathing 
mode was found to enhance electronic pairing [6] [23]. Fi- 
nally, experimental data on nanopowders [24] also indicate 
that one should not be overly pessimistic about cluster cal- 
culations. YBa2Cu306+a; powders consisting of 0.7 nm thick 
grains with a 11 nm radius show the same Tc as the bulk, 
indicating that going to the thermodynamic limit could not 
be essential to understand pairing. Cluster calculations can 
be relevant and physically insightful concerning the interplay 
of various degrees of freedom on pair structure and formation. 
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